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APPROXIMATE METHOD OF INTEGRATION OF LAMINAR 
BOUNDARY LAYER IN INCOMPRESSIBLE FLUID * 

By L. G. Loitsianskii 

Among all existing methods of the approximate integration of 
the differential equations of the laminar "boundary layer, the most 
■vriLdely used is the method based on the application of the momentum 
equation (reference l) . The accxiracy of this method depends on 
the more or less successful choice of a one-parameter family of 
velocity profiles. Thus, for example, the polynomial of the fourth 
degree proposed by Pohlhausen (reference l) does not give velocity 
distributions closely agreeing ■vri.th actual values in the neighbor- 
hood of the separation point, so that in the computations a strong 
retardation of the separation is obtained as compared with experi- 
mental results (reference 2) . The more-accurate methods employed 
in recent times (references 2 to 4) assume as a single -parameter 
family of profiles the exact solutions of some special class of 
flows with given simple velocity distributions on the edge of the 
boundary layer (single term raised to a power, linear f-unction) . 

The transition to the more complicated two- and more -parameter 
families of profiles would require, in addition to the momentum, 
equation, the employment of other possible equations (for example, 
the equations of energy (reference 5) and others (reference 6)). 

A greater accuracy might also then be expected for relatively simple 
velocity profiles that satisfy only the fundamental boundary con- 
ditions on the s\irface of the body and on the edge of the boundary 
layer. This second approach, however, as far as is known, has not. 
been considered except for very simple solution for the case of 
axial flow past a plate (reference 7). 

In the present paper, a solution is given of the problem of 
the plane laminar boundary layer in an incompressible gas; the 
method is based on the use of a system of equations of successive 
moments (including that of zero moment, the momentum equation) of 
the equation of the boundary layer. Such statement of the problem 


"Priblizhennyi Metod Integrirovania Uravnenii Laminarnogo Pogra- 
nichnogo Sloia v Neszhimaemom Gaze.” Prikladnaya Matematika i Mek- 
hanika, USSR, Vol. 13, no. 5, Oct. 1949, p. 513-525. 
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leads to a complex system of equations, which, however, is easily 
solved for simple supplementary assumptions . The solution obtained 
is given in closed form, by very simple formulas and is no-less 
accurate than the previously mentioned complicated solutions that 
are based on the use of special classes of accurate solutions Of 
the boundary -layer equations . 


1. Derivation of Fundamental System of Successive Moments of 
Boundary -Layer Equation. The well-kno’/'?n equations of the stationary 
plane laminar boundary layer in the absence of compressibility 
have the form 
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} ( 1 . 1 ) 
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vrhere u(x,y) and v(x,y) are the projections of the velocity at 
a section of the boundary layer on the ajcial and transverse axes of 
coordinates x and y, U(x) is a given longitudinal velocity on 
the outer boundary U' = dU/dx, and i? is the kinematic coeffi- 
cient of viscosity. \Vhen the equation of continuity is applied, 
the first of equations (l.l) m.ay be given the more convenient form 


L(u,v) = [u(U-uj] + [v(U-u)] + U'(U-u) 

ox . -* oy 


d^(U-u) ^ 


( 1 . 2 ) 

'The left side of equation (1.2) is multiplied by y^‘ and 
integrated mth respect to y from zero to infinity in the case 
of an asymptotically infinite layer or from zero to the outer 
limit of the layer y = 6(x) for the assijmption of a layer of 
finite thickness. In either case, the following expression is 
obtained: 
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It is assumed in this equation and in what follows that, in 
view of the very rapid approach of the velocity difference U - u 
to zero as y — all integrals with the infinite upper limit have 
a finite value. 


For k = 0, 




dx 


0 


u(U-u) 




where the magnitude 
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w 



( 1 .^) 


(1.5) 


represents the friction stress on the surface of the hody. 

Equation (1.4), the well-known impulse or momentum equation, 
is readily transformed into its usual form 


where 



For k = 1, a new equation of the 'first moment' is obtained 
from equation (1.3) 
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and, in general, for k > 2, the equations of successively increas- 
ing moments are obtained 
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y*^u(U-u) dy - k 
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y^"^(U-u) dy + U' y^(U-u) dy 

0 \Jo 


p »,6 


= k(k-l)u 


y^"^(U-u) dy 


(1.9) 
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In all these equations, the transverse velocity v(x,y) is 
assiimed expressed in terms of the axial u(x,y) from the equation 
of continuity. 


It is now assumed that the fetraily of functions 
u = (x, y; ^1, ^2^ • • • ^k^ 


( 1 . 10 ) 


satisfies the boundary conditions of the problem with k param- 
eters ^ 2 ^, . . . , which are functions of x, such that the 

k successive moments of equation (1.3) 


p“-,5 
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j'^L(u°, v°) dy 


( 1 . 11 ) 


iDecome zero. On the assumption that it is permissihle to pass to 
the limit k — it would then he possible to state that the function 

u(x,y) = lim u° [x, y> A^(x), >^2(^)» • * * » ^k^^^ 
k-»® 

( 1 . 12 ) 

\rith parameters ^2(^)^ . . • ^ ^k(x) satisfying the infinite 

system of equations 
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yj 


0 


v°) dy = 0 


(k = 0, 1, 2, . . . ) 


or, what is equivalent, systems (1.4), (1.8), eind (1.9) will be an 
exact solution of the fundamental system (l.l) for the assumed 
boundary conditions. 


For this solution, it is merely necessary to recall the kno^cn 
theorem that a continuous function, all successive derivatives of 
which are equal to zero, is identically equal to zero (reference 8). 


The question of the proof of the validity of this theorem is 
not considered in the case of an infinite interval or of an inter- 
val the boundaries of which are functions of a certain variable 
with respect to which the differentiation is effected. A certain 
construction, not based it is true on a. rigorous proof, of the 
solution of the problem ’.-rill be employed with the aid of the suc- 
cessive equations of the moments of the basic boundary-layer 
equation . 


2. Choice of Parameters of Family of Velocity Profiles at Sec- 
tions of Boundary Layer. Special Form of Equations of Moments. 

As is seen from the previously discussed considerations, the funda- 
mental difficulty lies in the choice of a feimily of velocity pro- 
files (l.lO) and the determination of the parameters of the 

family. One of the simplest methods of the solution of this prob- 
lem is indicated herein. 

In the converging part of the boundary layer, the velocity 
profiles at various sections of the layer are known to be almost 
sirriilar; the velocity profile is deformed mainly in the diffuser 
part of the boundary layer do^mstream of the point of minimum pres- 
stire. ' The deformation of the profile consists of the appearance 
of a point of flexure that arises near the surface of the body and 
then moves away from it as the separation point is approached. 

The presence of this deformation of the profile near the sur- 
face should greatly affect the magnitude t ^ proportional to the 
normal derivative of the velocity on the surface of the body; it 
will therefore diminish to zero as the point of separation is 
approached. The deformation of the profile will have a smaller 
effect on such integral magnitudes as 5* and 5** a.nd very 
little effect on magnitudes that contain under the integral 
sign fxmcticns that rapidly decrease as the surface of the body 
is approached. 
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For the parameters characterizing the effect of the deforma- 
tion of the velocity profile, it is natural to assume those magni- 
tudes that depend relatively strongly on the deformation of the 
velocity profile. With regard to the magnitudes that vary little 
with the deformation of the velocity profile, however, it is natiiral 
to assume that they do not depend on the chosen parameters. 

For the fundamental parameters determining a change in the 
shape of the velocity profiles, which may he called form parameters, 
the nondimens ional combination of the magnitudes 6* and &** 

'irtll be employed with the given functions U(x) and U'(x) and 
physical constants, namely, the parameters 



For the computation of the remaining magnitudes in the equa- 
tion of moments according to the assumption, the velocity profile 
will be assumed in a section of the boundary layer in a form that 
does not depend on the parameters f, and H: 

§ = ( 2 . 2 ) 

This assvimption permits, as will be subsequently seen, obtain- 
ing on the basis of very simple computations a sufficiently accurate 
solution of the boundary- layer equations for arbitrary distribution 
of the velocity on the edge of the layer. The transformation of 
equation (1.6) will now be considered. 

If the parameter ^ is introduced, then by equation (2.1), 


d6** U’S** 

dx ^ U 


(2+H) 
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It Is not difficult to obtain finally 


--IP' 


For the transformation of the left side of equation (1.8). 
the first integral can he written hy equation (2,2) (t^ = y/S**) 


yu(U-u) dy = ^19(1^^) dri = 

0 


2^**2 


where the magnit’jde equal to 


T|(?(lJp) dq 


represents a constant computed hy the given function 9 (ti). 


In order to compute the following integral, the transverse 
velocity V is first expressed hy the formula 

Pa a / \ 

V = - ^ dy = - — U tp dq 

Sx Sx V / 


= - U'5** 


<pdT) - U 


cpdT) - U 5** 9 — 
dx 


or, when it is noted that 


dx dx 




y d8** 


1 d5** 


.2 dx “ " ^ 5** dx 


the following expression is obtained: 


V = - U'5** 


9 dq - U 
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There is thus obtained 




v(U-u) dy = u2s*» 
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v(U-u) dy = 5** - H-UU' 5** ^ (2.7) 

^ dx ^ 


where Hg and denote the constants 

rn 




H2 = 
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cp diijd-cp) dT) 


.0 



q3dii)(l^) dTi 




( 2 . 8 ) 


Finally, the last integral in equation (1.8) is transformed into 

OC 

y(U-u) dy = H^U 5**2 (2.9) 

where the constant is equal to 


H4 = I Tl(l-q)) dq 


( 2 . 10 ) 
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By substituting the integrals obtained in equation (1.8), 


— - HpU^5** + E,UU*8**2 + H.UU'5**^ = UU 

dx ^ dx ^ 


( 2 . 11 ) 


or by replacing 6**^/u = f/U' by equation (2.1) and carrying out 
the transformation, 

{h - H)i = -J^y (2-12) 


When the nev constants are introduced, 

1 


a = 


b = 


H1-2H2 

= 1 - 1=2 




> (2.13) 




the equation of the first moment is reduced to the form 


M . £1 (a-bf ) + 2 : f. 
dx U ' U’ 


(2.14) 


The third equation is obtained from the system (1.9) by setting 

k = 2. 


i. I y2u(U-u)dy - 2 / yv(U-u)dy -i- U’ / y2(U-u)dy = 2u / (U-u)dy 


dx 


(2.15) 


There is obtained, as before. 
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where the constant He is equal to 

n=o 


% = 


Jo 


Ti^cp(l-cp)dTi 


Further, by analogy with equation (2,7) 
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where 


Poo 


Hg = 


Jo 


n \n 9 - I (pdq) (l^)dq 
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H7 = 




0 V J 


cpdq (l-cp)<iTi 
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) (2.19) 
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to 


The last integral on the left side of equation (2.15) is equal 



^ y^(U-u)dy = HqU ^Hq = J Ti2(l-<p)dq/ (2.20) 


The integral in eqiiation (2.15) on the right reduces to the 
unknown parameter H 


(U-u)dy = U 



;) 


dy = = U S**H (2.21) 


By substituting the expressions obtained for the integrals in 
the second-moment equation (2,15), there is obtained, after simple 
transformations , 

i (3H5-2H6) j| = - (is+lv^Haj f] + J (SHs-ZHg) 2^ f 


( 2 . 22 ) 
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The system of three equations (2,3), (2.14), and (2,22) has 
thus been established for determining the three xmknown msignitudes f, 
^ , and H, The solution of this system is now considered, 

3, Determination of the Constants Approximate Formulas 

for Parameters f, and H, For the determination of the numeri- 
cal values of the constants 3i, H2, , , ,, H9, the form of the 
function (il) must be known. The simplest velocity profile in 
the theory of the asymptotic boundary layer is the velocity profile 
in the sections of the boundary layer of the flow past a plate. 

The function cp(q) for this case can easily be determined from 
the generally known t able o f values of the velocity ratio u/U as 
a function of ? = y^U/i;x/2 . 


Superfluous computations may be avoided by noting that the con- 
scants to be computed are connected with one another by certain 
simple relations. 


First of all, from equations (2.3) and (2.14), 


(3.1) 


By setting f = C, there is obtained a = 2 ^q, where 

denotes the magnitude C computed for the plate (U* = 0, f = O). 
From t’ne definition of C ^.nd from the known relations for the 
plate , 



0.4408 


(3.2) 

Further, by comparing with one another the magnitudes Hi, H2, 



(3.4) 
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0, that is, the ratio b*/b** 
2.59 

It is then easy to obtain the value of b by equations (2,13), 
(3.3), and (3.4). 

b = = a (2 Hi+%-H 2+%-H2+ iSg^) 

= a (4 Hi- 2H2 + |Hq^) = 4 + | = 5.48 « 5.5 (3.5) 

When df/dx is eliminated from equations (2.22) and (2.14), 

1 3Hc:-2Hc r "1 

H = (H5+H7 + ^Hq) f + [1 - (2H1+H3+H4) fj 

= (H5+H7 + Ihq) f + I (3H5-2H6) (1 " 1 (3.6) 

By setting f = 0, 

I (3H5-2He) . ^ S!5.89 (3.7) 

The only magnitude that must be computed again from the table 
of values 9(q) is the magnitude Hg + Hy + B.q/Z, Numerical inte- 
gration gives 


where Hq is the value of H for f = 
for a plate is equal, as is known, to 

Ho = ^ = ^♦'^21 ^ 
^ b** 0.664 


H5 + Hy + I Hq = 24.73 
after vhich there is immediately obtained 


(3.8) 


H = 2.59 - 7.55 f 


(3.9) 



MCA TM 1293 


13 


Substituting this expression for H in equation (3,1) gives 
5 = 0.22 + 1.85 f - 7.55 (3.10) 

Finally, integrating the simple linear equation (2.14) gives 


f = 



nx 


u^"^(e )de = 




u^*^U)d4 


(3.11) 


Equations (3.9), (3.10) and (3.11) give the required solution. 

The simple, approximate solution Just obtained is now compared 
with the actual values. The almost complete agreement of the val- 
ues of f obtained with the first approximation (which is practic- 
ally the only one that is applied) of the preceding works (refer- 
ences 2 and 3) will be noted. The closed-form relation between 
and f likewise differs little from the corresponding tabulated • 
functions in the references cited. 

For comparison, the curves ^(f) and H(f) obtained accord- 
ing to the formxQas of reference 2 and by the formulas (3.10) and 

(3.9) are shown in figure 1. The results obtained will also be 
centred with the formulas of Wright and Bailey (reference 9). An 
approximate method of computation of the laminar boundary layer is 
proposed therein in which the equation of momentum (1.6) is employed 
with T.^ and 5** substituted by the formiOas for the flow past 

a plate. By expressing the results of Wright and Bailey in the 
parameters of the present report, the analogs of equations (3.9), 

(3.10) , emd (3.11) are obtained. 

H = 2.59 

5 = 0.22 + 4.09 f 

f = 0.44 2^ 

U 

It is easily seen that this formula for f corresponds to 
equation (3.11) for b * 1, The straight lines for ^ and H 
shown dotted in figure 1 indicate the considerable deviation of 
the formulas of Wright and Bailey from more accurate formiilas 
presented herein. 


\ (3.12) 
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For confirmation, the particular case of the laminar boundary 
layer corresponding to the so-called single-slope velocity distri- 
bution at the outer boundary of the layer U = 1 - x will be con- 
sidered. This case has been theoretically solved and an exact 
solution in a tabulated form (reference 10 ) is available. The 
results of the recomputation of these accurate solutions in the 
form assvuned by the parameters are given in figure 2. Also shown 
for comparison are the corresponding curves obtained by the pro- 
posed approximate method and by the method of Wright and Bailey. 

4. Possible Methods of Rendering the Foregoing Solution More 
Accurate . The method described in the preceding sections was based 
on the assumption of a slight dependence of on the form param- 

eters f, and H. This assumption may be eliminated and the 
method rendered more accurate, although it thereby becomes con- 
siderably more complicated. 

In order to discuss the possible generalizations of the method, 
the complete system of equations, for example, for the three- 
parsimeter case is written out* that is, a three -parameter family 
of velocity profiles is assumed in place of equation (2.2). 


S = f, £, H) (n = (4.1) 


By substituting this velocity profile in the system of the 
three equations of successive moments (1.6), (1.8) and (2.15), 
there is obtained a system of three ordinary nonlinear differential 
equations that determine the magnitudes of the parameters f, 
and H: 


2: JL M 

2 U» dx 



+ Hf 




(4.2) 


% - - 
2 


He 


Ko + 


= ^ [l . (2Ei.E3.H4)f] 4 2 : (Hi- 1 Hj) f 



f “ 

dx \ 

^ Sh / dx 

J fu 

1 TT \ 


(4.3) 
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(SHs-SHg) + 1 

f 

Sf / 


f ^ 

dx 

L 1 SH 5 I 

ix. tlH 

u' r 

, 4, „ 1 „ \ 

1 

\ 2 

dx 


J +- EqJ f 



I fr (3H5-2H6) f (4.4) 


In vhich, in addition to the previous notations, the following 
definitions are chosen: 



(l-(p)dT) 


(l-qp)dTi 




"A 


>(4.5) 
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It is noted that, in the system of equations (4.2), (4.3), 
and (4.4), and % are not constant magnitudes , as pre- 
viously, hut known functions of the form parameters f, C, and Hj 
the form of these functions depends on the chosen feimily of 
profiles (4.1). 

The equations (2.3), (2.14), and (2.22) earlier employed evi- 
dently represent a particular case of the system (4.2), (4.3), and 
(4,4) on the assumption that the family of velocity profiles at 
the different sections of the boundary layer has the form of equa- 
tion (2.2); in other words, these profiles are similar to one 
another. All values of are of course then equal to zero and 

is constant. 

The proposed method may be rendered considerably more accurate 
by assuming, for example, the single-parameter feunily of velocity 
profiles 

^ = CP(ti; f) (4.6) 


Then 









and the system of equations (4.2), (4.3) and (4.4) is transformed 
as follows: 


1 JJ ^ /g, 1 uu*' \ 

2 U’ dx ^ ‘ 2 u,2y 


f + Hf 





(2%+H5+H4)f] 


U2 

U' 



f 



- (hs+H,.! HgjfJ + i ^ (3H5-2H6)f 


(4.7) 


(4.8) 


(4.9) 



NACA TM 1293 


17 


Equation (4.8) can be given the form 

^ 1 - (2E]^+H3+H4)f ^ ^ ^1 ~ 2 ^2 ^ 

. ^ %*i ^2 + (Ki+d%/^f)f Hi - i Hg + (Ki+aHi/af)f 

Ct O' 

(4.10) 

which represents a generalization of equation (2.12) where equa- 
tion (4.10) approximates equation (2.12) because of the small change 
in with change in the parameter f and the smallness of the 

magnitude (KQ^+SHj^/Sf )f in comparison with - Hg/2, This gen- 
eralization permits obtaining the integral of equation (4.10) by 
introducing a correction to the solution of equation (2.12). 

By dividing both sides of equation (4.9) by the corresponding 
sides of eqixation (4.8) and thus eliminating df/dx, there is 
obtained 

H = (H5+H7+I H8)f + 


t 


2H, 




Kii + 



Hi 






[1 - (2Hi+H3+H4)f] 





- i (3H5-2He)Jf 

(4.11) 


By similar considerations on the smallness of the magnitudes 
(K4+I/2 ^H5/^)f in comparison with (3H5~2H0)/4 and of 

(K^+SEj/Sf )f in comparison with - H2/2 and on the slight 
variability of it may be concluded that the value of H 

determined by equation (4,11) is an improvement in the accuracy of 
the approximate value of H according to equation (3.6). 


It may be remarked that in this more accurate approximation 
there is no longer that universal relation between the pareuneters H 
and f, independent of the form of the fxmction U(x), character- 
izing the given particular problem. The presence in equation (4.11) 
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of a second term vith the factor UIP '/U' shows that in the more 
accxirate approximation the magnitude H in a given section of the 
layer depends not only on the value of the form parameter f in 
this section, as was the case in the rougher approximation of equa- 
tion (3.6) or (3.9), hut also on the value of the magnitude UU''/U'^ 
in the section considered, that is, on the values of the func- 
tion U(x) and its first two derivatives. It is readily observed 
that the second term on the right side of equation (4.11) will 
give a small correction to the solution (3.6) for relatively small 
values of the magnitude UU''/U'^. 

The same considerations hold for the expression for ^ , which 
may be obtained by substituting df/dx from equation (4.10) and 
H from equation (4.11) into equation (4.7): 




1 - (2%+H3+H4)f 


Hi - r H, 




[2 4 



+ 



+ 


2f + (H5+H7+ I Hg)f^- 



As is seen, in this new approximation, in contrast to the pre- 
ceding one, there is no universal relation between ^ and f . The 
presence of a term with the factor UU'*/U'^ makes the magnitude I 
depend not only on the value of the parameter f but also on the 
form of the function U(x) and its first two derivatives in the 
given section of the boundary layer. 


It is of interest to remark that in this approximation the 
position of the point of separation of the boundeiry layer, that is, 
the value of x = Xg for which ^ is equal to zero, will no 
longer be determined by some universal value of the form param- 
eter fg, but in each Individual case the value of x = Xg must 

be determined for which the right side of equation (4.12) becomes 
zero. 
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By assuming a particular form of a family of velocity profiles 
(4.5), employing, for example, the sets of velocity profiles 
applied in the previous investigations (references 2 to 4), the 
values of the functions and Kj^ eire determined} the form 

parameters f, and H, that is, the thickness of the momentum 
loss the friction stress and the displacement thick- 

ness 6* may then he found \vLthout difficulty. The solution 
of equation (4.10) and the determination of H and ^ hy equa- 
tions (4.11) and (4.12) offers no particular difficulty. Further 
improvement in the accuracy requiring the solution of a system of 
the type of equations (4.2), (4.3) and (4.4) is hardly of practical 
interest . 

In the previous discussion, the scheme of the asymptotically 
infinite boundary layer was used, but similar equations may be 
obtained adso for the case where the bo\mdary layer is assumed to 
be of finite thickness. 

The method here proposed may evidently also be applied to the 
case of the thermal boundary layer. The characteristic feature of 
the method for the cases of both the dynamic and the thermsil bound- 
ary layer lies in the fact that the friction stress and the quantity 
of heat given off by a unit area of the body sure expressed in inte- 
gral form and not in terms of the derivatives of functions that 
represent the approximate velocity and temperature distributions 
in the sections of the boundary layer. 


Translated by S. Reiss 
National Advisory Committee 
for Aeronautics . 


REFERENCES 

1. Loitsianskii, L. G.: Aerodynamics of the Boundary Layer. 

GTTI, 1941, pp. 170, 187. 

2. Loitsianskii, L. G.: Approximate Method for Calculating the 

Laminar Boundary Layer on the Airfoil. Comptes Rendus 
(Doklady) de I’Acad. des Sci. de L'URSS, vol. XXXV, no. 8, 
1942, pp. 227-232. 

3: Kochin, N. E., and Loitsianskii, L. G.: An Approximate Method 

of Computation of the Boundary Layer. Doklady AN SSSR, 

T. XXXVI, No. 9, 1942. 



20 


NACA TM 1293 


4. Melnikov, A. P.: On Certain Problems in the Theory of a Wing in 

a Wonideal Medium. Doctoral dissertation, L. Voenno- 
vozdushnaia inzhenernaia akademia, 1942. 

5. Leibenson, L. S.:. Energetic Form of the Integral Condition in 

the Theory of the Boundary Layer. Rep. No. 240, CAHI, 1935. 

6. Kochin, N. E., Kibel, I. A., and Roze, N. V.: Theoretical 

Hydrodynamics, pt. II. GTTI, 3d ed., 1948, p. 450. 

7. Sutton, W. G. L.: An Approximate Solution of the Boundary Layer 

Equations for a Flat Plate. Phil. Mag. and Jour. Sci., 
vol. 23, ser. 7, 1937, pp. 1146-1152. 

8. Carslaw, H., and Jaeger, J.; Operational Methods in Applied 

Mathematics . 1948 . 

9. Wri^t, E. A., and Bailey, G. W.: Laminar Frictional Resistance 

with Pressure Gradient. Jour. Aero. Sci., vol. 6, no. 12, 

Oct. 1939, pp. 485-488. 

10 . Howarth , L . : On the Solution of the Laminar Boundary Layer 

Equations. Proc. Roy. Soc. (London), vol. 164, no. A919, 

Feb. 1938, pp. 547-579. 





